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PACS 31.14.V - Electron correlation calculations for atoms, ions and molecules 
PACS 32 . 80 . -t - Photoionization and excitation 
PACS 34 . 80 . Dp - Atomic excitation and ionization 

Abstract. - The ionization dynamics of a two-electron atom in an attosecond XUV-infrared 
pump-probe experiment is simulated by solving the time-dependent two-electron Schrodinger 
equation. A dramatic change of the double ionization (DI) yield with variation of the pump- 
probe delay is reported and the governing role of electron-electron correlations is shown. The 
results allow for a direct control of the DI yield and of the relative strength of double and single 
ionization. 
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sfl Introduction. — With the emergence of coherent, ul- 
, ^ trashort laser pulses in the extreme ultraviolet (XUV) 
^regime the time-resolved investigation of ultrafast pro- 
^fcesses inside atoms became experimentally within reach, 
q often accounted for as 'attophysics' fil. Typically, 
— la higher-harmonics generated (HHG) sub-femtosecond 
XUV pump pulse triggers the electronic motion which is 
O then probed by an infrared (IR) or optical pulse with ad- 
(f-)justable delay. Owing to the weak intensity of the XUV 
l/~) pulse, usually strong probing pulses are needed, which 
Cn are nowadays routinely controllable and reproducible [2j . 
. These new techniques allow for the time-resolved explo- 
it") ration of fundamental sub- femtosecond dynamics in the 
^^ electronic properties of atoms [3||5], molecules [6] and 
^— ; atoms on surfaces [7] . 

• • Often the complicated treatment of multi-electron ef- 
^fects is simplified utilizing a single-active electron (SAE) 
description, together with the sudden approximations for 
S-Jthe XUV pulse. While this seems to yield satisfactory 
^ agreement with a certain class of recent experiments [8U9] , 
the validity range of the SAE remains open. As we will 
show, in many pump-probe scenarios (in particular, with 
intense probe pulses) two-electron effects are crucial and a 
full description of both, electron-electron (e-e) correlations 
and the two laser pulses, is necessary. 

In this Letter we address two important multi-electron 
processes in atoms: (i) the laser induced shake-up effect 
during XUV photon absorption and (ii) the rescattering 
mechanism. In (i) the XUV-photon removes one elec- 
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tron, causing a rapid change of the binding potential of 
the remaining electron(s). During their rearrangement, 
the XUV-photon energy is shared between the outgoing 
and the remaining electron(s) which is a consequence of 
e-e correlations. On the other hand, (ii) is a continuum 
effect and plays an important role for double ionization 
(DI) in the strong field regime. This mechanism has been 
investigated in great detail for helium, 



e.g. 



10 , and the 



importance of e-e correlations as well as the responsible 
mechanisms for non-sequential (NS) DI have been sub- 
ject of many studies, where finally experiments based on 
the COLTRIMS (cold-target recoil-ion momentum spec- 
troscopy 11 12]) technique favoured rescattering with im- 
pact ionization [l3}fl5] . However, still today not every 
facet of this process is understood 16 . Most importantly, 



the time dependence of nonsequential DI and, in partic- 
ular, its manifestation in a time-resolved pump-probe ex- 
periment, have remained unexplored. In this Letter, we 
study these questions. We demonstrate that mechanisms 
(i) and (ii) are of key importance for DI and that they are 
intimately coupled via electron energy transfer. Moreover, 
double ionization can be controlled to a high degree via 
mechanism (hrj 

Finally, we predict, that time-resolved observation of 
such processes is experimentally within reach. 



1 The controllability of mechanism (ii) through carrier-envelope 
phase modulation of a single pulse has be en d emonstrated in |17| and 
within a two-color few-cycle situation in [18| . However, we focus on 
control by delay and intensity in an XUV-IR two-pulse experiment. 
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Method. — The simplest system where the above 
mentioned effect are expected to occur is a two electron 
atom. Let us consider two electrons in a binding poten- 
tial Vbi and a time-dependent perturbing field V cx . Their 
motion is described by the two-particle time-dependent 
Schrddinger equation (TDSE), which reads 



ift*(r 1 ,r a ,t)=<|-i(^ 1 +^,) + 



(1) 



J2[V hi (r t ) + V CK ( ri ,t)} + w(\ ri ~ r 2 \)\ y{ ri ,r 2 ,t) 



i=l 



All quantities are given in atomic units (a.u.) throughout 
this work (h = m e = |e| = l/47reo = 1), unless stated 
explicitly. Equation (nl in its full dimensionality has been 



solved for various strong field situations, see e.g. 19 20 
and references therein. However, for an adequate treat- 
ment of the excitations considered in this work, the so- 
lution of the full two-electron problem is, yet, not man- 
ageable. Therefore, it is reasonable to consider a one- 
dimensional model atom, which has been successfully ap- 
plied to describe strong-field interaction with helium, e.g. 
21-23]. However, due to its long-range character the ID 
Coulomb binding potential still introduces high complex- 
ity in the description of the two-electron continua giving 
rise to a high computational effort (extremely large grids 
are required to simulate XUV-IR pump probe scenarios). 

In order to avoid these long-range binding effects, here, 
Vbi(Xj) is modeled by a ID potential well of width 2a 
and depth vo where Vbi = outside [—a, a]. We will 
show below (cf. Conclusions) that this model correctly 
reproduces key features of the two-color ionization dy- 
namics of ID helium. We use a depth of vq = 1.6 and 
a width of 2a = 5.6 which allows for a reasonable de- 
scription of a helium-like atom, including the ground state 
energy and a sufficient number of bound two-electron 
and single-electron (ion) eigenstates. A similar model 
has recently been used to compute the attosecond (as) 
electron pump-probe dynamics of surface atoms in solids 
[24p5 . In the spirit of the one-dimensional helium model, 
the electrons interact via a softened Coulomb potential, 
w(xi — x 2 ) = \j \J(x x — x 2 ) 2 + 1. With that, the ioniza- 
tion potentials for single ionization (SI) and DI evaluate 
to 4 1 ' = 0.92 and I ( p 2) = 2.41, cf. inset in Fig.[l] 

The two-color laser field is treated within the dipole 
approximation with Gaussian envelopes 



-x,{^ex P (-^ 
-E;° r exp ( -—2 ) cos(w ir t) 



cos[uj x (t — t)] 



*cx\Xij t) 



Throughout, we use a 240 as [a x = 10) pump pulse with 



r(i) 



(2) 



a photon energy of I p < hw x — 1.99 (54 eV) < Ip and 
intensity of 8.8- 10 13 W/cm 2 . These parameters are chosen 
such that significant SI is observed but multi-photon (MP) 
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Fig. 1: (color online) Kinetic energy distribution of an XUV ex- 
cited photoelectron from TDSE, (black) solid lines, and TDHF 
calculations, (red) dashed lines. The inset depicts the (exact) 
two-electron and the corresponding single-electron (ion) energy 
spectrum. The photon energy fiw x , schematically indicated by 
the (red) long arrow, is in between the first and second ioniza- 



r(l) 



,(2) 



tion threshold Ip and I p . The dominant shake-up process 
with energy transfer AE to the second excited ion state is in- 
dicated by the (gray) dashed arrow. 



absorption is suppressed. The probing pulse is given by a 
strong 900 nm few-cycle laser pulse, cf. Fig. [2j The delay 
time r of the XUV pulse is defined with respect to the 
maximum IR intensity and varied in temporal steps of the 
order of <j x . 

We solve Eq. (fTl) in coordinate representation within 
a grid-based finite-difference method [26] for a series of 
delay times r. Due to high kinetic energies occurring 
in our simulations, as a result of XUV photon absorp- 
tion and subsequent IR acceleration, large spatial grids 
are needed to avoid reflections at the boundaries. Typ- 
ically, we use a box size of at least —500 < x. L < 500 
with a minimum of 8192 grid points for each particle. 
To achieve proper time resolution (r-dependence) we have 
performed about 50 runs for two different IR pulses, see 
below. As initial state, we choose the spin-singlet ground 
state which possesses a symmetric spatial wave function, 
V(x 1 ,x 2 ,t) = y(x 2 ,x 1 ,t). 

The observables of interest are the SI and DI yield 
defined as the probability to find, respectively, one and 
two electrons outside a certain (sufficiently large) distance 
from the binding potential, see e. g. [27, 28 . Due to the 
absence of long-range binding and delocalized states in our 
model the present procedure is particularly well suited to 
compute SI and DI and to discriminate between both. Ki- 
netic energy spectra of photoelectrons coming from singly- 
ionized systems are obtained by projecting the final two- 
electron wave function onto plane waves (single-particle 
states), integrating out the bound part of the two-body 
wave function. 
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Fig. 2: (color online) Single and double ionization yield versus 
delay time r obtained from solutions of the TDSE and within 
TDHF, respectively, for an IR probing laser intensity of J = 
8.8 • 10 13 W/cm 2 , [pulse (i)]. The DI graphs are normalized to 
zero for large delays (sequential DI). Furthermore, the TDHF 
curve is scaled by a factor 10. The temporal laser intensity 
profile, oc |i5(t)| 2 , is sketched by the (red) dotted line in the 
lower figure. The maximum impact energy is reached by an 
electron created at tg and rescattering at t*. 



Results. — Let us first consider an IR intensity of 
8.8 • 10 13 W/cm 2 [pulse (i)]. The corresponding SI and 
DI yields for different delays r are shown in Fig. [2] We 
first note that SI directly follows the instantaneous inten- 
sity of the IR pulse at time of XUV excitation. This is 
not surprising, since the additional electrical field from 
the IR pulse increases the resulting total intensity of the 
combined IR and XUV pulse and the XUV photon energy 
YvjJx is sufficient to ionize the system independently of the 
IR field strength. The DI, however, shows a completely 
different r-dependence. While for r > it follows the IR 
intensity, for negative delays it strongly departs from the 
latter. In contrast to SI the curve is not symmetric with 
respect to r = 0, and even shows a significant maximum 
at the leading edge of the IR pulse, around r « —40. Since 

(2) 

hu x < Ip ' , the DI is truly caused by non-linear processes, 
either by (at least) two XUV photon absorption or through 
a combined action of both pulses. Since the two-photon 
absorption is a sequential effect, the DI should in that case 
follow the form of the SI yield, which is obviously not the 
case, cf. Fig. [2j Further, as the IR intensity is too low for 
DI through strong field tunneling or MP ionization, only 
the combination of both pulses can be responsible for DI. 

XUV induced shake up. — We first consider the 
case of XUV-only excitation. Two different methods are 
employed: the full solution of the two-particle TDSE, 
Eq. ([I]), and the solution within the time-dependent 



Hartree-Fock (TDHF) approximation 29 neglecting, by 
definition, all e-e correlation effects. For recent compar- 
isons of time-dependent mean-field calculations with the 



exact treatment of the TDSE in the context of strong-field 
interactions addressing limitations and improvements, see 
e.g. 



30 33 



The kinetic energy spectrum of photoelectrons coming 
from singly-ionized systems is shown in Fig. IT] Two 
distinct peaks separated by the photon energy fuv x can 
be identified in both spectra, TDHF and TDSE, at the 
expected positions of ftw x — Ip ~ 1, for absorption of a 
single photon, and 2hu> x — Ip « 3, for two-photon absorp- 
tion, respectively. Ground state e-e correlations become 
apparent in a small general shift between TDHF and 
TDSE spectra (£™se = -2.41 and £™hf = -2.35). A 
remarkable feature is the presence of a third pronounced 
peak at an energy of -Ekin ~ 0.28 in the TDSE calcula- 
tions, corresponding to slow photoelectrons. The missing 
kinetic energy, AE = -E s hakcup = 0.83, can be directly 
associated with the energy difference from the ion ground 
state to the second excited state. Therefore, and from 
the fact that this contribution is absent in TDHF, cf. 
Fig.[T] we conclude that this peak is direct evidence of the 
population of shake up states by the remaining electron 
in the iorQ see inset of Fig. [1] 



Semiclassical rescattering model (SRM). To 

shed more light into the underlying physics, we follow a 
classical analysis of the electron trajectories in the contin- 
uum 34 . In many cases, a (semi-) classical picture turns 



out to describe a variety of strong- field processes. We men- 



tion here the discussion of above-threshold ionization 35 



and NSDI |36 37 incorporating Coulomb effects 38 and 
using few-cycle pulses 39 40l. In this Letter we restrict 



ourselves to an adapted, simple version of the common 
three-step model: Consider an electron "created" at time 
to in the continuum with momentum po [to coincides with 
the delay time r] which follows from energy conservation: 



p = ±[2(hujx — Ip — -Eshakeup)] 1 ' 2 - Here we took into 
account a possible kinetic energy loss due to shake up of 
the second electron. 

Within this model, the electron which was excited by 
the XUV pulse with momentum po, acquires in the IR 
field with vector potential A(t) the momentum 



p(t, t ) =p Q + -A(t) =p Q - [ dt E(t) 
c Jt„ 

leading to the classical trajectory [x(to) = 0], 



(2) 



x(t,t ) =Po(t-t ) 



dt / diE(t) . (3) 

U) J U) 



The condition x{t r ,to) = 0, for t r > to, gives the 



2 In fact, a careful inspection of the DI in Fig. |2l reveals, that this 
shake-up state can be depopulated via IR probing, noticeable from 
the fact that the yield for r < —100 is larger (state depopulated by 
IR pulse) than for r > 100. This is in analogy to the experiment by 
Uiberacker et al. 3 where shake-up state population is probed in a 
similar way. 
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Fig. 3: (color online) Classical trajectories of an electron in the 
continuum created at time to during the IR probing field (ii), cf. 
Eq. pi). Dashed (blue) lines correspond to escaping electrons, 
(red and grey) solid lines to returning electrons. The probing 
IR field is indicated by the (black) dashed-dotted line. 



time of rescattering t r and the impact energy at t r : 

ESRM{tr\to) = P(tr,to) /2. 

Solutions of Eq. (I3J) for different times to, initial momen- 
tum po = V^Ekhi extracted from the kinetic energy distri- 
bution in Fig.JTland an IR intensity of 3.5- 10 14 W/cm 2 are 
shown in Fig. [3] Trajectories for returning electrons (solid 
lines) are allowed during the ascending cycle of the IR 
pulse, whereas electrons created at different times (dashed 
lines) escape. 

The analysis of Eq. pi reveals that, for IR pulse (i) , only 
electrons with po > and initial kinetic energy p^/2 < 0.4 
[cf. Fig. Ill can be driven back to the ion by the IR field 
thereby obtaining a kinetic energy EsRM(tr) sufficient for 
impact ionization. The energy Esrm is shown in Fig. [4] 
versus to an d exhibits a maximum around r = tg = —40. 



an electron "created" 
i.e. 



around tX will be accelerated 



In fact, 

until t*, i.e. for the entire half cycle of positive IR field 
strength, cf. Eq. (pi) and Fig. [2] Interestingly, the maxi- 
mum position of Bsrm(t) is very close to that of the max- 
imum of the DI yield observed in the TDSE simulations, 
cf. Fig. HI clearly supporting the rescattering mechanism. 
Figure HJ also shows that the value of .Esrm is sufficient for 
impact ionization from the shake up stata^l E^ n . 

Note that the dominant fraction of electrons (with en- 
ergies around 1, cf. Fig. [T]) escapes without rescatter- 
ing. Hence, the energy loss due to shake up of the second 
electron is crucial for impact ionization. To further ver- 
ify the rescattering mechanism we have performed a se- 
ries of TDHF calculations for the whole range of delays 
r. Within TDHF both, SI and DI, follow the IR inten- 
sity profile, cf. Fig. [2j in contrast to full TDSE solutions. 
Thus, e-e correlation effects are responsible for the ob- 
served r-dependence of the DI yield. 



3 The peak of -Esrm refers to electrons with an initial energy 
Pq/2 Ri 0.28 (vertical dashed line in Fig. lk, whereas electrons with a 
smaller initial energy will return with a larger value Esrm exceeding 
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Fig. 4: (color online) Impact energy (full lines) of the XUV- 
excited electron created at to = T for the IR intensities I = 8.8- 
10 13 W/cm 2 [pulse (i)] and I = 3.5- 10 14 W/cm 2 [pulse (ii)]. For 
comparison, the double ionization yield from TDSE simulations 
(cf. Figs. [2] and [5| is shown by dashed lines. Horizontal 
gray lines mark the possible binding energies of the remaining 
electron (eigenstates of ion), i. e. the minimal kinetic energy of 
the returning electron necessary for impact double ionization. 



Controlling double ionization. Now the ques- 

tion arises whether it is possible to further increase the 
DI yield relative to SI. From the above analysis it is ob- 
vious that DI could be enhanced if the dominant contri- 
bution of XUV photo-excited electrons, i. e. electrons 
with Skin (to) ~ 1 in Fig. Ill would be able to impact ion- 
ize the remaining electron. In fact, this is easily achieved 
by increasing the IR intensity, as we demonstrate below 
for / = 3.5 • 10 14 W/cm 2 [pulse (ii)]. As before, SI fol- 
lows the IR intensity profile, and DI shows a peak around 
r = —40 which, however, is now dramatically enhanced, 
cf. Fig. [5] With the present intensity increase by a factor 
of four, an enhancement of DI by more than two orders of 
magnitude is achieved, whereas SI is not affected signifi- 
cantly. The comparison with the classical impact energy, 
cf . Fig. |4j shows again an excellent agreement between the 
maximum positions of -Esrm and the DI yield. Also, the 
peak height of Esrm now by far exceeds even the ground 
state energy of the electron in the ion, thus now DI oc- 
curs from all ionic bound states. The complete absence of 
this peak in TDHF calculations (which predict a DI yield 
which is three to four orders of magnitude too small) , cf . 
Fig. [5j proves again that e-e correlations are the origin of 
this effect. 

Conclusions. — In conclusion, we have demonstrated 
that in XUV-IR pump-probe scenarios DI proceeds via 
a combination of electron shake-up and rescattering, at 
low IR intensity, or via rescattering at high IR intensity. 
By properly choosing the delay between the XUV and IR 
pulse and intensity of the IR pulse the DI yield can be var- 
ied within two to three orders of magnitude, reaching up 
to several percent of the SI yield. A simple physical expla- 
nation has been given which is straightforwardly extended 
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Fig. 5: (color online) Same as Fig. [21 but for an IR intensity 
of I - 3.5 ■ 10 14 W/cm 2 , pulse (ii). Additionally, the results 
for a ID helium model are shown for comparison in the upper 
figure. 



to other laser pulse shapes. We verified that our predic- 
tions for a ID model atom can be directly applied to small 
atoms. This is supported by a series of two-electron TDSE 
simulations for a ID helium atom for the parameters used 
in this letter. The results are included in Fig. [5] Obvi- 
ously, our main observation, the non-monotonic delay de- 
pendence of DI, is correctly captured by the model atom. 
This indicates that the correlation-dominated DI mecha- 
nism is present also in case of a Coulomb potential. On the 
other hand, there are quantitative differences (note the 11- 
fold increase of the peak height) which can be attributed 
to differences in the excitation spectra. A more detailed 
analysis of XUV-IR excitation of ID helium requires very 
large computational effort and is subject of ongoing work 
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Since the laser parameters used in our calculations 
are typical for current experimental conditions, we expect 
that time-resolved measurements of this process should be 
possible and allow to shed more light into one of the most 
intriguing manifestations of correlated electron dynamics 
in nature. 
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